Dynamics of dark solitons in elongated Bose-Einstein condensates 



A. E. Muryshev\ G. V. Shlyapnikov^'^'^^ Ertmer"^ K. Sengstock"^^, and M. Lewenstcin^ 
(1) Russian Research Center Kurchatov Institute, Kurchatov Square, 123182 Moscow, Russia 
(2) FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands 
(4) Laboratoire Kastler Brossel* , 24 rue Lhomond, F-75231 Paris Cedex 05, France 
(3) Institut fur Quantenoptik, Universitdt Hannover, D-30167 Hannover, Germany 
(4) Institut fiir Laserphysik, Universitdt Hamburg, 20355 Hamburg, Germany 
(5) Institut fiir Theoretische Physik, Universitdt Hannover, D-30167 Hannover, Germany 

We find two types of moving dark soliton textures in elongated Bose-Einstein condensates: non- 
stationary kinks and proper dark solitons. Tiie former iiave a curved notcii region and rapidly decay 
by emitting phonons and/or proper dark solitons. The proper moving solitons are characterized by 
a flat notch region and we obtain the diagram of their dynamical stability. At finite temperatures 
the dynamically stable solitons decay due to the thermodynamic instability. We develop a theory 
of their dissipative dynamics and explain experimental data. 
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Recently, several spectacular experiments have demon- 
strated the creation of vortices [Q-^ and dark solitons 
iQ-^ in Bose-Einstein condensates (BEG) of trapped 
alkaU-atom gases. These experiments open an unprece- 
dented possibility to study the dissipative dynamics of 
such macroscopically excited Bose-condensed states. The 
literature concerning the dynamical and thermodynamic 
stability of vortices is quite extensive (see, e.g., for 
review). As the vortex has a topological charge (circu- 
lation), the dynamically stable (single-charged) vortices 
can decay only when reaching the border of the conden- 
sate. At finite temperatures in non-rotating traps, the 
motion of the vortex core towards the border is induced 
by the interaction of the vortex with thermal excitations 
and is rather slow. The lifetime of vortices in trapped 
condensates jsj^ is thus relatively long and extends to a 
few seconds in the experiments |^|^. 

Dark solitons have a density dip and a phase slip in 
one direction and, as well as vortices, they are particular 
solutions of the Gross-Pitaevskii (GP) equation. Exten- 
sive studies of dark solitons in nonlinear optics [0 have 
expounded their transverse dynamical instability in 3D 
geometries, leading to the undulation of the soliton plane 
and decay into vortex-antivortex pairs and phonon waves. 
This scenario is similar to that observed at NIST Q for 
solitons in almost spherical trapped BEG's. Recently, the 
decay of solitons into vortex rings was observed at JILA 
Q. The transverse instability of dark solitons can be 
suppressed by a strong radial confinement of the soliton 
motion in elongated traps jl^. Dynamically stable soli- 
tons in such traps are not, however, thermodynamically 
stable, and their dissipative dynamics is expected to be 
fundamentally different from that of vortices. 

In contrast to vortices, the soliton has no topologi- 
cal charge and can decay without reaching the border 
of the condensate. Dark solitons behave as objects with 
a negative mass. The scattering of thermal excitations 
from the soliton decreases its energy, and the soliton ac- 



celerates towards the speed of sound, gradually looses its 
contrast, and ultimatelydisappears. This mechanism has 
been proposed in Ref . [IJ] , and the lifetime of the soliton 
has been obtained in terms of the reflection coefficient 
of the excitations. However, the theory of Ref. per- 
tains to the ID case, where the GP equation is integrable 
and the reflection coefficient is strictly zero within the 
Bogolyubov approach. Thus, one expects very long life- 
times of solitons in this limit. On the other hand, in 3D 
elongated traps the GP equation is no longer integrable 
and the scattering of thermal excitations from the soliton 
should be efficient. The absence of topological charge and 
integrability should then lead to a much faster dissipative 
dynamics of solitons than that of vortices. The Hannover 
results ||] indeed suggest that moving solitons generated 
in a cigar-shaped trap are dynamically stable, but their 
contrast decreases to zero within typically ~ 15ms as a 
result of the thermodynamic instability. 

In this Letter we study the dynamics of moving dark 
solitons in 3D elongated Bose-Einstein condensates and 
present three important results: i) We find that using 
the " phase imprinting method" |5|jl^ one can generate at 
least two kinds of soliton textures: non- stationary kinks 
and proper dark solitons. The former have a notch re- 
gion that moves with radially non-uniform velocity and 
undergoes bending similar to that observed at NIST [Q. 
These textures are dynamically unstable and decay via 
the emission of phonons and/or proper dark solitons. The 
proper solitons are characterized by a flat notch region 
and propagate without changing their shape; ii) We de- 
rive the diagram of dynamical stability for proper soli- 
tons; iii) We solve the problem of reflection of excitations 
from the soliton and analyze its decay due to thermody- 
namic instability. The dissipative dynamics exhibits an 
interplay between the extent of non-integrability and the 
absence of topological charge, and the soliton lifetime 
ranges from milliseconds for Hannover-type 3D solitons 
to more than seconds in quasi ID geometries. 
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We consider a condensate with repulsive interaction 
(the scattering length a > 0). The condensate wave func- 
tion can be written as ^'(r, t) exp {—int), where fi is the 
chemical potential. In an infinitely long cylindrical har- 
monic trap this function satisfies the GP equation 

^^^4-|^A + f .y+,|*(r,Ol^v,(r,.). (1) 

Here ujp is the frequency of the radial (p) confinement, 
g = 47r?i^a/m, and m is the atom mass. The wave func- 
tion of the ground-state condensate minimizes the corre- 
sponding energy functional and is the solution of Eq. (|^) 
with zero Ihs. Macroscopically excited Bose-condensed 
states (solitons, vortices etc.) are described by any other 
stationary, or time dependent solution of Eq.(|l|). Obvi- 
ously, they are thermodynamically unstable as they do 
not correspond to the minimum of the energy functional. 

A stationary, or solitary-wave macroscopically excited 
BEG state can also be dynamically unstable with regard 
to elementary excitations around it. The unstable excita- 
tion modes are characterized by complex eigenfrequencies 
and grow exponentially in time, which indicates that the 
BEG state will evolve far from the initial shape. 

Strictly speaking, dark solitons are solutions of the ID 
GP equation in free space. They are characterized by a 
local density minimum (notch) moving with a constant 
velocity v, and by a phase gradient of the wave function 
at the position of the minimum. In an otherwise uniform 
condensate of density uq, the dark soliton state is de- 
scribed by the wave function (see ||l^ and refs. therein) 

vl/(z,t) = y/n^{cos0 - isin6'tanh[sin6'(z - wt)/^o]) , (2) 

where cos^ = v/cs, and Cs = \J nQg/m is the speed of 
sound. The quantities and —9 are the phases of the 
soliton state at {z — vt) — )■ — oo and (z — vt) oo, re- 
spectively. The width L of the notch (soliton plane) is of 
the order of the correlation length Iq — h/mcs- Approxi- 
mate soliton solutions can be also found in ID harmonic 
traps [Ti|JT5[ | . Solitons which have non-zero velocity in 
the center of the trap oscillate along the trap axis [g5| . 

In 3D harmonic traps the solutions of the GP equation, 
describing standing dark solitons {v = Q and d'^ /dt = 0), 
have been found in |Tl|,|I^. For infinitely long cylindri- 
cal condensates these solutions follow from Eq.(||) and 
can decay due to the transverse dynamical instability. 
The stability criterion requires a strong radial confine- 
ment providing a non-Thomas-Fermi (TF) regime with 
the radial size of the condensate r L ^ Iq. 

The existence of moving soliton-like textures in 3D 
elongated condensates is confirmed by the experiments 
and simulations |^], but no analytical solution has been 
found so far. In the TF regime (/x 3> JitOp), one can use 
Eq.(H) with p-dependent uq and Iq. Then, the absence of 
the radial flux of particles at an infinite axial separation 
from the notch requires the phase 9 to be independent 
of the radial coordinate. This means that the notch ve- 
locity V depends on p and is proportional to the local 



velocity of sound Cs{p)- Hence, the central regions of 
the notch move faster than the borders, and an initially 
flat notch region starts bending in the course of motion. 
Our simulations for TF condensates show that such non- 
stationary kinks can be created by phase imprinting with 
/9-independent optical potential. They are dynamically 
unstable and decay on a time scale of the order of ujp^- 
The notch surface bends more and more, whereas the 
notch velocity increases and the depth decreases. Ulti- 
mately, the non-stationary kink decays emitting excita- 
tion waves and (for p,/fujjp < 10) a proper dark soliton: 
The latter moves with p-independent velocity and does 
not change its shape characterized by a flat notch region. 
In non-TF condensates {p^fiujp) we generated the proper 
solitons directly by simulating the phase imprinting. 

In the absence of dissipation, the velocity u of a proper 
soliton remains constant in an infinitely long cylindrical 
condensate. The wave function depends on p and x = 
z — vt. In order to find this wave function, we write it 
in the form \I'(p, x) — 'ip{p, x)f{x), where the functions 
'ip{p,x) and f{x) satisfy the equations 

ihd^P/dt = {~{h^/2m){A+ 2(V,///)V,)+ mLjy/2 

+9\fm'-m\)}^, (3) 

ihdf/dt - -ih^/2m)A^f + (pif) - p)f, (4) 

The quantity p is a functional of / and has to be found 
self-consistently from Eqs. (^) and (^). We will select 
the function ip{p,x) such that at infinite x it becomes 
the wave function of the ground-state condensate, ipQ^p). 
Hence, for \x\ oo we have |/| 1 and p ^ p. 

We consider the limiting case where the axial size L 
of the soliton notch greatly exceeds the radial size r of 
the condensate. Then the radial distribution of parti- 
cles is close to no{p) = "00 (p) for the ground-state con- 
densate. The quantity p{f) is close to p and can be 
expressed as p{f) = p + (|/P — l)gdp/dg + Sp, where 
(5/2 is a correction of higher order in r/ L. The quantity 
gdp/dg = mc^, where Cs is nothing else than the ve- 
locity of axially propagating sound waves in the ground- 
state condensate. In the quasi ID regime, where the 
interparticle interaction at maximum condensate density 
nomd ^ ^i-^p, we have an almost Gaussian density pro- 
file no(p) — riom exp (— p^/Zp). The radial size r ^ Ip — 
(?i/mwp)^/^, and the small parameter of the expansion 
for p is [rjV)^ ~ n^ragl'hijJp- To first order in n^^raglfiujp 
we obtain gdp/dg — nomg/2. This gives the velocity Cg 
which is by a factor of -^2 sma ller than t he speed of sound 
at maximum density: Cg — ^/nomS/Sm. For radially TF 
condensates the chemical potential p cx ^Jg and we arrive 
at the same expression for Cg- This result for TF elon- 
gated condensates has been obtained in and found in 
the MIT experiment jlj]. The condition r <^ L requires 
fast TF solitons for which the density dip is small and 
the function |/(a;)| is close to 1. 

Omitting the higher order correction 5p, Eq.(^) for the 
function f{x) becomes an ordinary ID GP equation 
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iWf/dt - -{h^/2m)^J + mc2[(|/|2 _ 1)/. (5) 

The terms in the rhs of Eq.(||) are related to the axial 
kinetic energy and to the mean-field interparticle inter- 
action in the presence of the radially inhomogeneous den- 
sity profile of the Bose-condensed state. The dark soliton 
solution f{x) is then given by the rhs of Eq.(||), where 
Cs is replaced by Cs, by To = and no by unity. 

The soliton-like wave function of the condensate can be 
written as ^'(p, a;) — ipo{p)f{x). This clearly shows that 
the maximum soliton velocity (at which the soliton dis- 
appears) is equal to Cs, i.e. is by y/2 smaller than the 
sound velocity at maximum condensate density. 

In the quasilD regime, to first order in riomS/^Wp the 
function ip{p^x) = ipoip) + S'ip{p,x), where a small term 
(5-0 is real. The second order correction Sfl is then equal to 
5(1/1' - l)?9'd't^/dg\ We find g^d^p/dg' = -jnomg, 
with 7 = 37iom5'lii(4/3)/4?iWp <JC 1. Thus, we obtain 
A(./) = Ai+mc2[(|/|2-l)-7(|/|2_i)2], and Eq.(|) becomes 



the reference frame moving together with the soliton, 
£ — hkv = e' — fik'v. For small k the refiection coeffi- 
cient increases as fc^ oc e^. The coefficient reaches its 
maximum at fc ~ 1^^ yj\ — v'^ /c^ and decays exponen- 
tially for large k. The calculations leading to Eq.(^ are 
quite lengthy and will be published elsewhere. 

The reflection of excitations from the soliton provides 
a momentum transfer from the thermal cloud to the soli- 
ton. Hence, there is a friction force acting on the soliton. 
The momentum transfer per unit time is given by 



h{k - k')R(k)v{k)N{e - hkv)dk/2TT, (8) 



dt 2m, dx^ 



The quasilD dark solitons are dynamically stable, and 
we now discuss their thermodynamic instability in the 
presence of a thermal cloud. With 7 = 0, Eq.(^) trans- 
forms to Eq.(||) which is integrable. Hence, the dark 
solitons described by this equation are transparent for 
thermal excitations. There is no energy and momentum 
exchange between the soliton and the thermal cloud, and 
the thermodynamic instability does not manifest itself. 
The dissipative dynamics of the solitons originates from 
the interaction between the radial and axial degrees of 
freedom, which in the quasi ID regime is described by the 
small term 7(|/P — 1)^ in the rhs of Eq.(^. This term 
only slightly modifies the wave function vE'(p, x), but it 
lifts the integrability of the equation and is responsible 
for the reflection of excitations from the soliton. 

For finding the refiection coefficient we have solved the 
Bogolyubov-de Gennes equations following from Eq.(||). 
For the phonon branch of the spectrum, where the ax- 
ial momentum of an excitation k <C Iq^ , the excitation 
wave functions u, v were found in the form of expansion 
in powers of fc^o and 7 around the fundamental modes of 
the Bogolyubov-de Gennes equations with 7 = 0. The 
w, V functions of these equations were obtained straight- 
forwardly for an arbitrary k and used for calculating the 
refiection coefficient R{k) from the Fermi golden rule. 
For klo <C 1 the obtained R{k) matches the one following 
from the method of fundamental modes. At any energy 
e and momentum k of the incident wave we have 



where N(e — hkv) are equilibrium occupation numbers 
for the excitations. 

The energy of the soliton can be written in the form 
H = {Mc^j3){l-v^ /clf/^ , with M = 2TTno^lor^m being 
the effective mass of the soliton {r — \p2lp in the quasilD 
regime). The soliton energy decreases with increasing 
V. For example, for w <C Cs we have H — Mcl/i — 
1) -7(1/1' - (6) Mv'^/2. The quantity = M/m = 2imQjQr'^ > 1 

is the number of particles that one has to remove from 
the condensate in order to create the soliton density dip. 
Thus, the dark soliton can be treated as a heavy classical 
particle-like object with a negative mass, and the friction 
force accelerates the soliton towards the velocity of sound 
(see ||l^]). The Hamiltonian equation dH/dp — v gives 
p = —Mv{l — v'^/c^Y^'^, and using Eq.(||) we obtain the 
time dependence of the soliton velocity. 

For T 3> nomg the time t at which the soliton contrast 
C — {1 — I'^/cg) decreases from the initial value Co to 
C{t), can be found from the relations 



F{C) - F{Co) = t/T, T = hN,/TRo, 



(9) 



R{k)= 



8TT"f{e — hkv)ll 



9?isinh{(7r(|fc|+|fc'|)?o/2Vl-u2/c2} 



\kk'\ 



v{k)u{ky 



(7) 



where v{k) = d{e — T%kv)dfik is the group velocity. The 
energy e' and momentum k' of the reffected wave are 
related to e and k by the energy conservation law in 



where F{C) is a universal function of the contrast, and 
Rq = O.O8472 is the maximum value of the reffection 
coefficient in the limit of w — > 0. The function F{C) 
was calculated numerically and can be approximated as 
F{C) « 0.471n{(l - C)/C} The quantity r can 

be regarded as a characteristic lifetime of the soliton. 
For example, if the contrast is initially equal to 30%, it 
decreases to 10% at a time t ~ 0.25r. 

The dissipative dynamics of quasilD solitons is gov- 
erned by the small extent of non- integrability of Eq.(|^) 
(7^1) and the time r can be very long. The physi- 
cal picture changes to the opposite one if the axial size 
L of the notch becomes comparable with the radial size 
r of the condensate. Then the non-integrability of the 
GP equation is essential and the absence of topologi- 
cal charge provides a fast dissipative dynamics. In this 
case, which corresponds to the border of dynamical sta- 
bility of the soliton, the developed analytical approach 
can be no longer used and we have found the proper soli- 
ton solutions numerically. To investigate the dynamical 
stability of moving proper solitons we have numerically 
solved a time-dependent equation for elementary exci- 
tations around the obtained soliton wave function. For 
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a given soliton velocity v, the axial size of the notch de- 
creases with increasing the ratio nomg/fi^p- Above a crit- 
ical value nQmg/fu^p — the transverse instability was 
manifesting itself in our calculations as a dramatic rise of 
excitation modes. In Fig.l we present the critical ratio 
as a function of the soliton velocity. For nQmg/fiuJp < 
the solitons arc dynamically stable. Note that £,c — 2.5 
for a standing soliton, which agrees with the earlier cal- 
culation For the solitons with high velocities the 
stability condition is more relaxed. 

For large nQj^g/TiLOp, i.e. in the TF limit, the soliton 
velocity v required for the transverse dynamical stabil- 
ity approaches the sound velocity Cg. If nQrng/fi^p > 10, 
then the transverse instability is suppressed only for v 
larger than the Landau critical velocity calculated for 
elongated condensates in Solitons with u > are 

characterized by the longitudinal instability related to 
the Cherenkov radiation of axially propagating excita- 
tions. We have found this instability numerically and es- 
tablished that it develops on a time scale longer than that 
of the transverse instability. Thus, for n^rng/^p > 10 
dark solitons are always dynamically unstable. 



In conclusion, we have investigated the dynamical sta- 
bility and dissipative dynamics of solitons in elongated 
EEC's, and explained the experimental data of Ref. 
For recently achieved quasilD EEC's, the parameter 
7^0.1 and our theory predicts the soliton lifetime larger 
than seconds. This opens prospects for studying dissipa- 
tive phenomena originating from the quantum character 
of the boson field omitted in the common GP approach. 
Temperature dependence of the soliton lifetime offers in- 
teresting possibilities of EEC thermometry. 

We acknowledge support from Deutsche Forschungs- 
gemeinschaft, the Dutch Foundations NWO and FOM, 
Alexander von Humboldt Stiftung, the Russian Founda- 
tion for Fundamental Research. We thank L. Carr, J. 
Dalibard, P. Fedichev, E. Phillips, and A. Sanpera for 
discussions. 

*LKE is an unite de recherche de I'Ecole Normale 
Superieure et de I'Universite Pierre et Marie Curie, as- 
sociee au CNRS. 




0.2 0.4 0.6 

v„/c. 



FIG. 1. The critical ratio nomg/hujp = versus the soliton 
velocity v (in units of Cj,). Solitons are stable below this curve. 

The condition L ^ r was fulfilled for solitons in the 
Hannover experiment For the number of atoms 

iV « 1.5 X 10^, and the trap frequencies uo ^ — 2tt x 14Hz 
and u)p — 2tt X A25i{'L, we calculate the critical tempera- 
ture Tc ~ 350nK, the maximum density ~ 4 x 10^^ 
cm~'^, and the chemical potential /x « 140 nK. This in- 
dicates that the solitons were in the TF regime, with 
[ijhLOp K, 7. The thermal fraction was about 10%, which 
corresponds to T « 0.5Tc and /i/T « 0.8. The soliton 
contrast was decreasing from approximately 30% to be- 
low the resolution limit (10%) at a time of 15 ms. The 
results in Fig. 1 indicate that the dark solitons of Ref. H] 
were dynamically stable. Ey using the numerical tech- 
nique which will be described elsewhere, we calculated 
the reflection coefficient of excitations from the soliton. 
The dependence of i? on /c and v is similar to that in 
the limit of r L. The maximum reflection coefficient 
for v^O is Ro w 0.7. Then from Eq.(|) we find r « 80 
ms and conclude that the soliton contrast decreases from 
30% to 10% at a time of 20 ms, in agreement with H. 



[8 

[9: 

[lo; 

[11 
[12; 

[is: 

[14 

[is: 

[16: 
[17: 



[is: 

[19 
[20 
[21 



M.R. Matthew et ai, Phys. Rev. Lett. 83, 2498 (1999). 
K.W. Madison et ai, Phys. Rev. Lett. 84, 806 (2000); F. 
Chevy et ai, ibid. 85, 2223 (2000). 

C. Raman et al, Phys. Rev. Lett. 87, 210402 (2001). 
J. Denschlag et ai, Science, 287, 97 (2000). 

S. Burger et al., Phys. Rev. Lett. 83, 5198 (1999). 
B.P. Anderson et al, Phys. Rev. Lett. 86, 2926 (200 1). 
A.L. Fetter and A.A. Svidzinsky, |cond-mat/0102003| ; R. 
J. Donnely, Quantized Vortices in Helium //(Cambridge 
University Press, Cambridge, 1991). 

P.O. Fedichev and G.V. Shiyapnikov, Phys. Rev. A60, 
R1779, 1999. 

A.A. Svidzinsky and A.L. Fetter, Phys. Rev. Lett. 84, 
5919 (2000). 

Y.S. Kivshar and B. Luther-Davies, Phys. Rep. 298, 81 
(1998). 

A.E. Muryshev et al, Phys. Rev. A. 60, R2665 (1999). 

P.O. Fedichev et al, Phys. Rev. A. 60, 3220 (1999). 

L. Dobrek et al, Phys. Rev. A60, R3381 (1999). 

R. Dum et al, Phys. Rev. Lett. 80, 2972 (1999). 

Th. Busch and J. Anglin, Phys. Rev. Lett. 84, 2298 

(2000). 

D. L. Feder et al, Phys. Rev. A62, 053606 (2000). 

E. Zaremba, Phys. Rev. A57, 518 (1998); S. Stringari, 
ibid. 58, 2385 (1998); G. Kavoulakis and C. Pethick, ibid. 
58, 1563 (1998). 

M.R. Andrews et al, Phys. Rev. Lett. 79, 553 (1997). 
This expression is not valid for u ~ 0, where the soliton 
motion starts in the diffusive regime, see [Q. 
P.O. Fedichev and G.V. Shiyapnikov, Phys. Rev. A63, 
045601 (2001). 

A.Gorlitz et ai, Phys. Rev. Lett. 87, 130402 (2001); F. 
Schreck et al, ibid. 87, 080403 (2001). 



4 



